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Abstract 

\ We present a unified derivation of Bohmian methods that serves as a common starting point for 

the derivative propagation method (DPM), Bohmian mechanics with complex action (BOMCA) 
and the zero- velocity complex action method (ZEVCA). The unified derivation begins with the 

iS 

ansatz ip = where the action, S, is taken to be complex and the quantum force is obtained by 
writing a hierarchy of equations of motion for the phase partial derivatives. We demonstrate how 



different choices of the trajectory velocity field yield different formulations such as DPM, BOMCA 



and ZEVCA. The new derivation is used for two purposes. First, it serves as a common basis for 

OO : 

comparing the role of the quantum force in the DPM and BOMCA formulations. Second, we use 

in 

the new derivation to show that superposing the contributions of real, crossing trajectories yields 
a nodal pattern essentially identical to that of the exact quantum wavefunction. The latter result 

r» : 

suggests a promising new approach to deal with the challenging problem of nodes in Bohmian 
mechanics. 

C3 1 PACS numbers: 
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I. INTRODUCTION 



The contrast between the inherent non-locality of quantum mechanics and the locality of 
classical mechanics has driven a decades long quest for a trajectory-based formulation of 
quantum theory that is exact. In the 1950's, David Bohm, building on earlier work by 
Madelungl] and de Broglie^, developed an exact formulation of quantum mechanics in 
which trajectories evolve in the presence of the usual Newtonian force plus an additional 
quantum forcej^, 4|. Bohm's formulation was originally developed as an interpretational 
tool, to recover a notion of causality in quantum mechanics. In 1999 Lopreore and Wyatt[5| 
demonstrated that the Bohmian formulation can also be used a numerical tool to do quantum 
calculations. This innovation has coincided well with the ongoing interest of the chemical 
physics community in finding effective numerical tools for performing multidimensional quan- 
tum calculations. The apparently local dynamics of the Bohmian trajectories suggest the 
possibility of computational advantages compared to fixed grid methods or direct-product 
basis set methods that scale exponentially with dimensionality. The Lopreore and Wyatt 
paper has motivated the development of a variety of new numerical approaches for imple- 
menting Bohmian mechanics j^, Q, [s[ 10, 11, 12]. Reference [l3|] gives an excellent account 



of the progress in Bohmian related formulations in recent years. 

The Bohmian formulation has two main limitations that currently limit its usefulness as 
a numerical tool. First, note that the non-locality of quantum mechanics does not disappear 
in the Bohmian formulation — it manifests itself in the quantum force term, which needs to 
be calculated to propagate the quantum trajectories. The quantum force is very unstable 
numerically when the wavefunction is oscillatory. This is related to the second limitation: 
the quantum force diverges at nodes of the wavefunction resulting in the break down of 
the Bohmian formulation in the vicinity of nodes. Indeed, most of the new contributions to 
Bohmian methodology have been aimed at overcoming these two limitations of the Bohmian 
formulation. 

Among the approaches that have been developed to deal with the calculation of the 
quantum force are the derivatiye propagation method (DPM)[l^] and Bohmian mechanics 



with complex action (BOMCA) la. Il6l|. The two approaches use a similar procedure writing 



a hierarchy of equations of motion for partial derivatives of the phase. The difference between 
the methods lies in the use of a complex action in BOMCA as opposed to a real amplitude 



2 



and phase in DPM. This difference has far-reaching consequences that we explore in this 
paper. 

In the first part of this publication we present a unified and compact derivation of 
Bohmian methods that serves as a common starting point for the DPM and BOMCA meth- 
ods, and the Zero- Velocity Complex Action (ZEVCA) method. The derivation is similar to 
the derivation of BOMCA but it leaves the choice of the trajectory velocity field undeter- 
mined. We demonstrate how different choices of the velocity field yield different formulations, 
such as DPM, BOMCA and others. The unified derivation allows for ready comparison of 
DPM and BOMCA, and the role that the quantum force plays in both formulations. In 
the second part of this paper we use the new derivation to show that superposing the con- 
tributions of real, crossing trajectories yields an interference pattern. This contrasts with 
a cardinal principle of the exact Bohmian formulation, that trajectories are not allowed 
to cross in configuration space. However, since we are dealing with an approximation to 
Bohmian mechanics, the no-crossing role does not apply. Superposing the contribution from 
real crossing trajectories represents a promising new avenue for dealing with the challenging 
problem of nodes in Bohmian mechanics. 

This paper is organized as follows. In section [Til we present the unified derivation. Section 
[TTTJis dedicated to several special cases of the derivation: ZEVCA (IIIIAj) . BOMCA (IIIIBj) 
and DPM flHI Cj) . In Section [TV] we show how a variation of the derivation leads to an accu- 
rate description of interference and nodal patterns. Section [V] is a summary and concluding 
remarks. 



II. UNIFIED DERIVATION OF BOHMIAN METHODS 



We start by inserting the ansatz 13, [l8|, [19] 

ip(x, t) = exp 

into the TDSE 



h 



S{x,t) 



ifttpt 



2m 



Ipxx + V(x,t)l/J, 



(2.1) 



(2.2) 



where S(x,t) is a complex function, h is Planck's constant divided by 2tt, m is the mass 
of the particle and V(x,t) is the potential energy function. The subscripts denote partial 
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derivatives. The result is a quantum complex Hamilton- Jacobi (HJ) equation [171. 1 la. Il 



s t + —S 2 X + v = —s xx , 

2m x 2m 



(2.3) 



where we recognize on the LHS a complex version of the classical HJ equation. On the RHS 
is an additional non-classical term that we refer to as the "complex quantum potential". 
This term is different from the quantum potential in the conventional Bohmian formulation. 
The time-independent version of ansatz (12.11) is familiar as the starting point of the WKB 



approximation [20 



21 



22j. Equation (12. 3ft was derived by Pauli[l3] as a first step in the 



WKB derivation. More recently, eq. (I2.3H was rediscovered by Leacock and Padgett 23], 
who also reverted to a time-independent formulation to calculate eigenvalues. The equation 
has also been used by several other authors}^. 25, 26] as an analytical tool, but not as a 
constructive method to solve the TDSE with trajectories. 

The classical HJ equation can be conveniently solved by integrating along trajectories 
that satisfy the classical equations of motion. From a mathematical point of view, these 
trajectories are the characteristics of the classical HJ equation. Here we use an analogous 
approach to solve the quantum complex HJ equation, by integrating along some family of 
trajectories. We define a family of trajectories by choosing a velocity field v(x, t), which can 
be done in an infinite number of ways. (As opposed to the classical case, the velocity field 
is not necessarily predetermined by the PDE we are trying to solve.) Solutions of 



- = ,(x,t) 



(2.4) 



determine trajectories, parametrized by their initial position x(0). The change of the phase 
along the trajectory is given by operating on S(x, t) with the Lagrangian time derivative 



with the result 



d _ d dx d d d 
dt dt dt dx dt dx 1 



—— = S t + vS x = -^-S xx Si — V + vS x 

dt 2m 2m x 



(2.5) 



(2.6) 



where we have used eq. (12.31) . The integration of eq. fl2.6f) requires S x [x(t),t] and S xx [x(t),t], 
i.e. the values of the first two spatial derivatives of S along the trajectory. Fortunately, it is 
possible to write equations describing the evolution of the spatial derivatives of S. Writing 

d n S 



S n [x(t),t] 



dx n 



(2.7) 



[x(t),t] 



and taking the n spatial derivative of eq. fl2.3p . we have 

(S t ) n + ^(SX + V n = ^S n+2 . (2.8) 

Inserting the result in the definition of the Lagrangian time derivative of S n gives 

dS 1 ih 

-7T = (St)n + vS n+1 = -—(Sf) n -V n + vS n+ i + —S n+2 . (2.9) 
at 2m 2m 

From this procedure we obtain an infinite set of coupled equations describing the evolution 

of S and its spatial derivatives along a trajectory, namely eqs. fl2.9p for n — 0, 1, 2, We 

note that evaluation of S n requires knowledge of S n+2 and S n +i (the latter through the term 
(Sf) n = 2^=0 \j)Sj+iSn-j+i)- The fact that evaluation of S = S requires knowledge of 
£>2 — S xx expresses the nonlocality of the Schrodinger equation. However from a formal 
point of view, integration of the infinite hierarchy of eqs. (12.91) along the trajectories defined 
by eq. fl2.4p can be regarded as a local method of solution of eq. (12.31) in the sense that infor- 
mation propagates along individual trajectories independently. This does not contradict the 
nonlocality of the Schrodinger equation, since not only is the value of S being propagated 
down a trajectory, but also all its spatial derivatives S\, S 2 , ■ ■ ■ are being propagated. 

A numerical approximation for solving eq. (12.31) can be obtained by truncating the set 
( 12. 9p at some n = N, by setting SV+i = Sn+2 = 0. We summarize the equations of motion 
of the approximation: 

^ = v[x(t),t], (2.10) 
^ = -— {Sl) n -V n + vS n+l + ^-S n+2] n = 0,l,...,iV, (2.11) 
Sn+i = Sn+2 = 0, (2-12) 

where we emphasize that the solution is given along an individual trajectory x(t) determined 
by the velocity field v (x, t). The initial conditions of eqs. (12. Ill) are given by 

<9 n ln[^(x,0)] 



S n [x(0),0] = -ih 



dx r - 



(2.13) 

x(0) 



where we have used the relation S(x, 0) = — ih\\i[%p{x, 0)] from ansatz (12. ip . The wavefunc- 
tion at time tf at final position x(tf) is given by 

^[x(t f ),t f ] = expUs [x(t f ),t f }\. (2.14) 



Equations fl2.10l) - fl2.12l) provide a common starting point for deriving several quantum tra- 
jectory methods such as DPM, BOMCA and ZEVCA as well as adaptive grid techniques. 
These methods differ from each other by the specific choice of the velocity field. For the 
sake of simplicity we present and compare the different methods for N = 2. For this case, 
eqs. (12.1 01) - (12. 121) yield four equations of motion 

£ = (2,5a) 

f-il-^+J^ < 2 - 15b » 

^1 = _^s 2 - V, + vS 2 , (2.15c) 
at m 

^ = -^-V 2 . (2.15d) 
at m 

The N = 2 case is unique for two reasons. First, it is the lowest order of truncation for which 
the equation of motion for the phase So (eq. fl2.lip for n = 0) includes a quantum potential 
term ^S 2 . Second, eqs. fl2TT5D for BOMCA, ZEVCA and DPM yield the exact solution for 
an initial Gaussian wavepacket propagating in a potential with up to quadratic terms. 



III. SPECIFIC CHOICES OF TRAJECTORY VELOCITY FIELDS 

A. Zero- Velocity Complex Action (ZEVCA) 

The simplest choice of the velocity field is 

v (x, t) = x(t) = x(0). (3.1) 

The resulting trajectories are straight lines, hence we refer to the resulting approximation 
as the Zero- Velocity Complex Action method[27] (ZEVCA). The ZEVCA formulation can 
be regarded as a hybrid between a grid method and a local semiclassical method. In ref- 



erence 



131 ] section 7.2, Wyatt considers the solution of the global hydrodynamic equations 



of quantum mechanics on fixed grid points (Eulerian grid) but dismisses its usefulness as a 



numerical tool. In reference 27J the ZEVCA formulation is shown to produce useful output 



from local propagation at a single grid point. In this paper we focus on the relation between 



BOMCA and DPM, hence we wil 
interested reader is referred to 27]] 



not elaborate further here on the ZEVCA method. The 
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B. Bohmian Mechanics with Complex Action (BOMCA) 



In the BOMCA method, the velocity of the trajectories is set as 

v\x(t),t] = SMt) ' t] . (3.2) 
m 

The rationale for this choice is evident if we recall that Si has units of momentum and that 
So can be associated with a quantum action field (eq.f l2.3l) ). Inserting eq.f l3.2j) in eqs.f l2.15l) 
yields the N = 2 equations of motion of BOMCA 

dx Si 

— = — , (3.3a) 
at m 

f = -K, ^ (3.3c, 

^ = -Vt- (3.3d) 
at m 

The BOMCA formulation was originally presented in terms of equations of motion for partial 
derivatives of the velocity instead of the phase [ljj. The transformation between the two 
formulations is straightforward by using relation (13. 2p . 

Equations (13.31) have several attractive features. 1) Equations (13.3a ) and (13.3b ) are New- 
ton's second law of motion in disguise, resulting in characteristics that are classical trajec- 
tories. In other words, the N = 2 BOMCA approximation does not incorporate a quantum 
force in computing the trajectories. Note, however, that generally the classical trajectories 
take on complex values of position and momentum. This results from the complex initial 
momentum that emerges from eq. fl2.13p for n = 1. 2) Recognizing the RHS of eq. (13.3b ) as 
a "quantum Lagrangian" , this equation has the familiar structure of the equation of motion 
for the classical action, with the addition of the quantum potential. The quantum potential 
is the only /i-dependent term, hence this term is entirely responsible for incorporating the 
quantum effects in the N = 2 BOMCA approximation. Note that for N = 3, a quan- 
tum force term appears in the trajectory equations (specifically eq. fl3.3b )). yielding complex 
quantum trajectories. 

Equations (I3.3a -d) appear in the context of other time- dependent semiclassical meth- 
ods that use complex classical trajectories. The first of these is the generalized Gaussian 



wavepacket dynamics (GGWPD) of Huber and Heller [28|, |29(. The two main advantages 



of the BOMCA formulation over GGWPD is that 1) the latter does not provide correction 
terms to the N = 2 approximation, and 2) the latter cannot be generalized to a non-Gaussian 
initial wavefunction. The second context in which eqs. (I3.3|) appear is in the first-order for- 



mulation of the complex trajectory method (CTM) of Boiron and Lombardi|30|, [32j. This 
formulation is a complex trajectory version of time-dependent WKB. Although the CTM 
method is identical to BOMCA for N = 2, it differs at higher orders of N. 

The propagation of complex trajectories requires the analytical continuation of the initial 
wavefunction to the complex plane, as well as a method to reconstruct the wavefunction on 
the real axis at the desired final time. Regarding the latter, in reference 15|] we describe 
an algorithm to calculate complex initial positions that end at a final time tf on the real 
axis = 0). Without going into the details of the algorithm, suffice it to say that it 

exploits the local analyticity of the mapping x(0) — ► x(tf). By inserting the phase So that 
corresponds to these trajectories into the original ansatz (12.11) . we obtain the wavefunction 
at a real position x(tf) at time tf. 



C. Derivative Propagation Method (DPM) 

We present here a simplified derivation of the DPM method that yields an e quiv alent but 
more compact set of equations than has appeared previously in the literature [13|. Setting 
the velocity field to be 



v[x(t),t] 



^{SMt),t]} 



m 



and inserting eq. (13.41) into eqs. fl2.15l) . yields the N = 2 DPM equations of motion: 

> 

[S(Si)] 2 



dx 

~dt 

dSo 



m 



dt 
dSi 

~dT 
dS 2 
~~dT 



2m 

m 



V + ^S 2 
2m 



2m 



-Vo 



si 

m 



(3.4) 

(3.5a) 
(3.5b) 
(3.5c) 
(3.5d) 



The trajectories in the DPM propagate on the real axis, removing the need to extrapolate 
the wavefunction to the real axis, as is required in BOMCA. But the classical structure that 
was evident in BOMCA is no longer present: comparing eqs. (13.31) and (13.50 . we see that a 
quantum force term appears in the equation of motion for the momentum S\ (eq. (l3.5b )). 
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and the equation of motion for the phase Sq (eg. (13. 5b )) has an "extra" quantum potential 



term ^' ;Sj ^ 



Equations (I3.5P are equivalent to those of the DPM as they appear in the literature {13I fl^| . 
as we now show. In the DPM, the wavefunction is represented by the ansatz 



ip(x, t) = exp 



c(x,t) + ^s(x,t) 



(3.6) 



where s(x, t) and c(x, t) are real functions (ansatz (I3.6P is essentially identical to the conven- 
tional Bohmian ansatz ip(x, t) = A(x, t) exp [|s(x, t)\ if one identifies A(x, t) = exp[c(x, t)}). 
Equating eq. (13.61) with eq. (l2.ip yields 

S(x,t) — s(x,t) — ihc(x,t). (3.7) 

Inserting eq. (13.71) into eqs. (l3.5p and dividing the results into their real and imaginary parts 
yields 

£ = sr" + 5S<* + * (3 - 8b) 

-1 = -Vy + — C1C2, 3.8c 

at m 

— = + — c 2 - V 2 , 3.8d 

at mm 

dc s 2 . . 

-dt = " W (3 ' 8e) 

dCl ClS2 (3.8f) 



dt m 
dc 2 2c 2 s 2 



(3.8g) 



dt m 

These equations are equivalent to eqs. (l3.5a -d). and are readily seen to be equal to eqs. (10.10) 
of reference [l3] by taking N = 2 and identifying x(t) = 51 . 

Note that the conventional DPM equations of motion also have a quantum correction 
to the trajectories at truncating order N = 2. This can be seen in the equation for s% 
(eq. (l3.8b )). as well as in the equations for c , c\ and c 2 (eqs. fl3.8e -g)) which have no classical 
counterparts. Even though the quantum force term in eq. (13.8b ) is proportional to ti 2 , a 
simple example demonstrates that this term is not small in any way compared with the 
classical force. Consider the harmonic oscillator potential V(x) = ^mu 2 x 2 , where uj is the 
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angular frequency. The time-evolution of the ground state is 



tf)(x, t) = exp 



2 + h 



(3.9) 



where a = ^ff and 70 = -f hi(^). We recall that eqs.flSJU) (and eqs.([33D) yield the exact 
solution for an initial Gaussian wavepacket in a harmonic potential. Comparing eq. ()3.6p 
and eq. (13.90 shows that s is position-me/epene/ent, hence, si = 0. As a results, according 
to eq.( 13.8h ) the eigen-trajectories are straight lines, a result familiar from the conventional 



Bohmian formulation. Since -J 1 = 0, we conclude from eq. (13.8b ) that the quantum force is 
equal in its magnitude to the classical force 

Vi[x(0)] = mw 2 i(0). 



h 2 

— cic 2 
m 



(3.10) 



The LHS of this equation is the quantum force. Not only is the quantum force not negligible, 
we see from the RHS of eq. fl3.10p that it depends linearly on the position. The appearance 
of a significant quantum force in this most "classical" example can be viewed as a result of 
working with real trajectories. In the BOMCA formulation with its complex trajectories, the 
quantum force vanishes in this case, which is one of the strongest motivations for studying 
BOMCA|15]. 



IV. INTERFERENCE WITH REAL TRAJECTORIES 



As mentioned in the Introduction, the nodal problem is currently the main obstacle to 
performing numerical calculations using Bohmian mechanics. In this section, we present 
preliminary results showing that an oscillatory wavefunction in close agreement with the 
quantum result can be obtained using eqs.f l2.15l) with real trajectories. Consider an initial 
Gaussian wavepacket 



ip(x, 0) = exp 
propagating in a Morse potential 



-a (x - x c 



+ tPc{x - x c ) + -70 



(4.1) 



V{x) = A{[1 -exp(-/3x)] 2 - 1} . 



(4.2) 



The parameters of the initial Gaussian are «o = 0.5, x c = 9.342, p c = 0, 70 = — ^ln(-^ a ) 
where we take m — h — 1 (all quantities are given in atomic units). As for the Morse 
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FIG. 1: An initial Gaussian wavepacket propagating in a Morse potential. The parameters of the 
system are given in the text. The final propagation time tf equals roughly half of the oscillation 
period of a classical particle located at the center of the initial Gaussian, x c . The exact wavefunction 
was calculated using the split operator method with a Fourier basis. 

parameters, A = 10.25 and (3 = 0.2209. The final propagation time, tf — 5.93, is chosen so as 
to produce a strongly oscillating pattern. Figure [1] depicts the potential and wavef unctions. 

In both DPM and BOMCA the velocity is taken to be a function of S\. in DPM, v = 
9fc ) while in BOMCA, v = We now present an alternative procedure; we take v to 
be completely independent of S% and set a classical equation of motion for v. 

dv _ Vxjxjt)] 

dt m 1 ' 1 

Equation (14. 3[) is then supplemented to the set of eqs. fl2.15l) . Furthermore, we extend the 
freedom in the choice of the velocity field to include the initial conditions for the velocity, 
which we take to be 

+(0 ),0] = ?»M3[ I (0)- I ,]; (4.4) 
m m 

Si[x(0),0] is obtained by inserting eq. fl4.ll) into eq. fl2.13l) and setting n — 1. Equation 
f)4.4p defines real initial conditions; taken together with eqs. (I4.3|) and (12. 15a ) this yields 
real classical trajectories (note that if we had taken the initial velocity to be complex, 
t>[x(0),0] = ; we would have obtained the BOMCA equations, eqs. fl3.3p ). The next 

11 



X 



FIG. 2: Real classical trajectories obtained by solving Newton's second law of motion (eq, (|4.3p ) 
for a Morse potential (eq. (|4.2p ) with initial conditions (|4.4|) . The parameters of the propagation 
are given in the text. The final propagation time tf is marked explicitly. The trajectories can 
be divided into two overlapping branches. The first branch (marked as solid lines) is the locus of 
trajectories that have reached their classical turning point and were reflected by the exponential 
barrier of the potential. The second branch (marked as dashed lines) is the locus of trajectories 
that by tf did not reach their classical turning point. 

step is to solve eqs. fl2.15p and (14.31) with initial conditions ( 12. 13[) and (14.41) . respectively, for 
a series of initial positions {x(O)} in the vicinity of x c . The wavefunction at time tf at final 
position Xf is given by eq. (l2.14[) . where S [x(tf),tf] is the solution of eq. (l2.15b ). 

In figl2]we plot the trajectories obtained. The trajectories can be divided into two over- 
lapping groups that we refer to as branches. The first (reflected) branch (marked as solid 
lines) is the locus of trajectories that have reached their classical turning point and were 
reflected by the exponential barrier of the potential. The second branch (marked as dashed 
lines) is the locus of trajectories that by tf did not reach their classical turning point. Thus, 
to an arbitrary final position x (x > —2.8) these correspond two initial positions and two 
associated trajectories. 

The question arises, if two trajectories and therefore two values of So(x,tf) correspond 
to each final position x, how should we determine the wavefunction ip(x,tf)7 The clear 
distinction between the two branches allows us to associate a wavefunction with each branch, 
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which we will call i/;i(x,tf) and if)2(x,tf). We apply the superposition principle to reconstruct 
the wavef unction at x, obtaining 

i>(x,tf) = if)!(x,tf) +tl> 2 (x,tf). (4.5) 

In figj3]we compare the exact wavefunction along with \ipi(x, tf)\, \tp2(x,tf)\ and \ip(x,tf)\ = 
\ijji(x,tf) + ip2{x,tf)\. For values of x larger than m —2.8 (the position of the maximum of 
\ip(x,tf)\), the superposition yields a surprisingly accurate approximation of the oscillating 
wavefunction: even though the wavefunctions ipi(x,tf) and ip2(x,tf) exhibit no oscillations 
whatsoever, their superposition yields strong oscillations and a node near the maximum. 
Note that ip 2 {x,tf) provides the "main" contribution to the final wavefunction; ipx(x,tf), 
which originates from the reflected branch, contributes to ip(x,tf) only where the wavefunc- 
tion oscillates. As the classical turning point is approached, the amplitudes arising from 
the different branches diverge and the superposition of contributions appears to have no 
physical significance. 

As noted above, the BOMCA equations for N = 1 are identical to those of Generalized 
Gaussian Wavepacket Dynamics (GGWPD). In both cases, the trajectories that are prop- 
agated obey classical equations of motion but are complex. In this section we described a 
modification of BOMCA in which the classical trajectories are taken to be real. It is in- 
teresting to speculate if there might be a connection with thawed Gaussian propagation, in 
which the equations of motion are the same as those for GGWPD, but the trajectories are 
real. Consistent with this conjecture is the observation that in thawed Gaussian propaga- 
tion there is no need for a root search, as is true in the real-trajectory version of BOMCA. 
However, the correspondence cannot be exact. In the real-trajectory version of BOMCA, in 
principle every point in coordinate space is propagated, whereas in thawed Gaussian prop- 
agation the initial wavefunction is decomposed into Gaussians and any decomposition that 
satisfies completeness is allowed. Moreover, in real-trajectory BOMCA the trajectories have 
no width and no functional form whereas in thawed Gaussian propagation there is always 
some residual signature of the Gaussian decomposition. For example, Gaussians that reach 
the turning point have part of their amplitude extending into the classically forbidden region, 
whereas in real-trajectory BOMCA such penetration into the classically forbidden region is 
absent. 
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FIG. 3: A comparison between the exact wavefunction with the contribution of the first branch 
\ipi(x,tf)\, the second branch \ip2(x,tj)\ and a superposition of both \ijj(x,tf)\ = \ipi(x,tf) + 
ip2(x,tf)\. For x > —2.8 (the position of the maximum of the exact wavefunction) the superposition 
procedure yields a surprisingly accurate approximation of the oscillating wavefunction: even though 
the wavefunctions \ipi(x,tf)\ and \ip2{x,tf)\ exhibit no oscillations whatsoever. The procedure is 
undefined in the classically forbidden region which the classical trajectories cannot access. The 
exact wavefunction was calculated using the split operator method with a Fourier basis. 

V. SUMMARY 



We have presented a unified and compact derivation of Bohmian methods that serves as a 
common starting point for several approximations. In particular, the new approach was used 
to derive and compare the DPM, BOMCA and ZEVCA methods. We focused on the role of 
the quantum force in the DPM and BOMCA methods, and showed that for N = 2, the lowest 
order of truncation for which a quantum potential is present, the BOMCA formulation is 
closer to the classical equations of motion than DPM, although the trajectories are complex. 
We also demonstrated that an interference pattern can be obtained by superposing the 
contributions from real classical trajectories in what is otherwise essentially a Bohmian 
formulation. 

Since nodes in quantum mechanics arise from interfering amplitudes, it is only natural 
to attempt to solve the nodal problem in Bohmian mechanics by applying the superposition 
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principle. This suggests decomposing the wavefunction into two nodeless parts and propagat- 



ing each part separately using trajectories 
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3l| . However, since nodeless wavefunctions do 



not generally remain nodeless, this approach generally requires a series of time-dependent de- 
compositions of the total wavefunction. Such decompositions are numerically inconvenient, 
largely arbitrary, and typically valid only for relatively short times. In the scheme we pre- 
sented, an oscillatory wavefunction is successfully decomposed for all times considered into 
two nodeless, non-oscillatory parts, ip\ and ife. As such, the scheme represents a promising 
new avenue for dealing with the challenging problem of nodes in Bohmian mechanics. We 
are currently exploring an analogous idea of applying the superposition principle to complex 
crossing trajectories in BOMCA 16J] and in time-dependent WKB32(; the results are very 
promising and will be published elsewhere. In addition, although we do not yet have final 
results, it seems that path integral techniques can provide a rigorous justification of the 
superposition principle in Bohmian methods and the limits of its applications. 
This work was supported by the Israel Science Foundation (576/04). 
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